Velocity relaxation in a strongly coupled plasma 
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Collisional relaxation of Coulomb systems is studied in the strongly coupled regime. We use 
an optical pump-probe approach to manipulate and monitor the dynamics of ions in an ultracold 
neutral plasma, which allows direct measurement of relaxation rates in a regime where common 
Landau- Spitzer theory breaks down. Numerical simulations confirm the experimental results and 
display non-Markovian dynamics at early times. 



More than half a century ago, Landau [T] and Spitzer 
[5] derived simple expressions for Coulomb collision 
rates that have become fundamental to modern plasma 
physics. Precise knowledge of collisional relaxation rates 
is essential for understanding plasmas of all varieties. It 
is fundamental to energy exchange in multi-species sys- 
tems [2] and determines transport properties, such as self- 
diffusion rates as well as thermal and electric conductiv- 
ities [3]. The underlying theory, however, breaks down 
in strongly coupled systems such as Jovian planet interi- 
ors |3] and dense-plasma experiments [S], which display 
strong correlations between particles. Here, we present 
the first direct measurement of thermalization rates in 
an unmagnetized, strongly coupled plasma. Exploiting 
the very low temperatures in ultracold neutral plasmas 
[5J[7], we realize strong coupling conditions at low enough 
densities to enable direct time-resolved measurements via 
optical manipulation and imaging. The observations are 
supported by numerical simulations that moreover high- 
light the importance of non-Markovian relaxation effects. 

In weakly interacting systems, that are either very hot 
and/or very dilute, relaxation is dominated by binary 
small- angle scattering events of distant particles. Conse- 
quently, a test charge traversing a single-species plasma 
of temperature T and density p, undergoes Brownian mo- 
tion with a corresponding damping coefficient [TJ [5] 
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where the factor TZ(v) derives from the so-called Rosen- 
bluth potential [5] and m denotes the mass of the test 
particle and the plasma charges. The term In A, known as 
the Coulomb logarithm, is determined by an upper cutoff 
for possible impact parameters that ensures convergence 
of the relaxation rate. In the original Landau-Spitzer 
derivation it is set equal to the Debye screening length, 
beyond which interactions are collectively screened by the 
surrounding plasma charges. 

Equation is applied to a wide range of plasmas, 
but it is only valid when spatial correlations in the sys- 
tem are weak. The degree of particle correlations can be 
characterized by the ratio of their average potential and 
thermal energy, as expressed by the Coulomb coupling 
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FIG. 1: (color online), (a) Schematics of the experimen- 
tal approach to probe ion relaxation in an ultracold plasma. 
Two counter-propagating, circularly polarized lasers, detuned 
by A PP /27T = -20 MHz from the 5s 2 Si /2 - 5p 2 Pi /2 transi- 
tion, optically pump population between the two ground-state 
magnetic sublevels (m = ±1/2) of ions in an ultracold stron- 
tium plasma. The corresponding level-scheme is shown in (b) , 
also indicating excited-state decay with the spontaneous emis- 
sion rate 7. The optical pumping produces skewed velocity 
distributions f±(v) for each of the ground states, which we 
probe via ion fluorescence induced by a circularly polarized 
light sheet, applied at a variable time t after optical pumping. 
A typical, simulated time evolution of the velocity distribu- 
tions, /+ (red) and /_ (blue), during the optical pumping 
stage (gray) and subsequent relaxation is shown in (c). 
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where a = (4irp/3)~3 is the average distance between 
the plasma charges. A plasma becomes strongly cou- 
pled when T > 1, i.e. when interactions start to dom- 
inate thermal motion. The Coulomb logarithm can be 
written in terms of the coupling parameter, A ~ T -3 / 2 , 
showing that In A turns negative for V > 1, and the 
Landau-Spitzer rate ([!]) becomes entirely nonsensical in 
the strongly coupled regime. 

Theoretical efforts to understand relaxation under 
strong-coupling conditions have largely focused on dense 
plasmas [9HI3]. as produced by intense-laser heating of 
solid-state samples [14] . An accurate description of re- 
laxation and transport processes is essential for the in- 
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terpretation of these experiments [TS] and, in particu- 
lar, for optimizing conditions for inertial confinement fu- 
sion |16j . Experimental probes of dense plasmas have 
advanced tremendously [171 I18j . and allow indirect in- 
ference of relaxation rates through theoretical modeling 
of other observables [THI [20]- However, direct and pre- 
cise measurement of relaxation rates remains an open 
challenge [21] . largely due to the fast dynamical time 
scales, complicated initial conditions and complex evo- 
lution at solid density. Ultracold neutral plasmas [5J [7] 
present an appealing platform for studying strongly cou- 
pled plasma physics under simple and well-controllable 
conditions. Because of their low densities, ultracold neu- 
tral plasmas evolve slowly enough that the dynamics of 
many plasma parameters can be measured directly [221 - 
125] with high temporal resolution. 

We create an ultracold neutral plasma by photoion- 
izing laser-cooled strontium atoms in a magneto-optical 
trap. Peak plasma density is varied from p ~ 10 9 — 
10 10 cm -3 by changing the delay between release of the 
trapped atoms and photoionization. By tuning the ion- 
ization lasers, we set the initial electron temperature to 
T e (0) = 105 K while the initial kinetic energy of the ions 
equals the thermal energy (~ lOmK) of the laser-cooled 
atoms. This yields weakly coupled electrons (T ~ 10~ 2 ) 
and would place the ions deep into the strongly cou- 
pled regime with V ~ 10 3 . However, photoionization 
produces initially uncorrelated ions, such that the sub- 
sequent development of ion correlations results in strong 
heating [3D] to T ~ IK during the first few 100 ns [5T1 - 
134] , and yields ionic Coulomb coupling parameters of or- 
der unity. Ion temperatures are determined by fitting the 
Doppler-broadened laser-induced- fluorescence spectrum, 
while the ion density is determined from absorption mea- 
surements of the plasma ions [35] . 

Due to the vast electron-ion temperature disparity, 
electronic screening of ion-ion interactions is considerably 
weaker than direct screening by the ions. For our condi- 
tions, the electronic Debye length is about 5 to 10 times 
larger than average particle distance a. The ionic com- 
ponent can thus be viewed as a classical one-component 
plasma, where the electrons provide a neutralizing back- 
ground and electron-ion collisions play a negligible role. 
The latter drive the relaxation dynamics of dense plas- 
mas [HHH]. Here, the time scale for electron-ion relax- 
ation greatly exceeds the duration of our measurement, 
which allows us to exclusively study ion-ion relaxation 
processes in the following. 

To this end, we exploit the degeneracy of the Sr+ 
5s 2 S 1( / 2 ground state, i.e. the availability of two distinct 
electronic spin states (see Fig. [lji). Initially, both spin 
states are equally populated, but we can manipulate their 
populations by applying two counter-propagating laser 
pulses with identical frequencies, detuned by A pp /27r = 
—20 MHz, opposite circular polarizations, and peak sat- 
uration parameters so = 6 (see Fig. [IJl). Due to the 
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FIG. 2: (color online), (a) Ion fluorescence spectra of a 
plasma with density of p = 10 9 cm -3 and temperature of 
T = 2.2K, corresponding to a Coulomb coupling param- 
eter of r = 1.2. The circles show the measured spectra 
from the m = +1/2 states at different times t after op- 
tical pumping. The spectra without optical pumping are 
shown by the triangles. They match the familiar Voigt profile 
arising from the underlying Maxwellian velocity distribution 
feq(v z ) = f+(v z ) + f-(v z ) of all plasma ions. The lines are fits 
of equation |3| to the experimental spectra, which yields the 
velocity distributions shown in (b) . Panel (c) shows the corre- 
sponding average velocity (dots), along with the result of our 
quantum-classical simulations. The squares are obtained by 
calculating fluorescence spectra from the simulated velocity 
distributions, which are then analyzed just as the experimen- 
tal spectra are. The line shows the average velocity obtained 
directly from the simulations. The good agreement between 
both approaches confirms the accuracy of our experimental 
procedure for measuring the ion velocity. 

Doppler effect, this optically pumps population between 
the two spin states in a velocity-selective manner. Ions 
with velocities around v z = A pp /fc along the wave vec- 
tor k of the a~ laser are transferred from m = +1/2 
to m = —1/2 and vice versa around — v z . This pro- 
duces skewed velocity distributions f±(v z ) for each spin 
state as shown in Fig. [ij. Note, that the state of the 
plasma and in particular the total velocity distribution 
fi v z) = f+(t>z) + f-( v z) remains undisturbed and pre- 
serves its Maxwellian form. Hence, the relaxation of f± is 
driven by an equilibrium plasma with a well defined tem- 
perature T. In this way, our experiments realize the orig- 
inal Landau-Spitzer construction of tagged test charges 
evolving in an equilibrium plasma background [2J. 

At an adjustable time after optical pumping, we ap- 
ply a third, less intense probe beam (so = 1.5) with a~ 
polarization and an adjustable detuning A pr and record 
laser-induced fluorescence spectra perpendicular to the 
beam propagation (Fig. [iji). We shape the probe beam 
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into an ellipsoid that only interacts with a narrow cen- 
tral sheet of the plasma, and, thus, produces an im- 
age of a two-dimensional cut through the plasma cloud 
as shown in Fig. [IJi. We analyze fluorescence from a 
small central area where the density is nearly constant 
and the hydrodynamic expansion velocity 29 is negligi- 
ble. The spectrum is given by a convolution of f+(v z ) 
and the Lorentzian profile of the probe transition, and 
can, thus, be used to determine the ionic velocity dis- 
tribution. Figure [2] shows a typical sequence of spectra 
obtained at different times after optical pumping. Ini- 
tially, one observes a considerable asymmetry with popu- 
lation enhancement around A pr = — A pp as expected. At 
later times, the spectrum approaches the familiar Voigt 
profile, indicating relaxation towards a Maxwellian ve- 
locity distribution within a few 100 ns, which is close 
to the inverse of the corresponding plasma frequency, 
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collisional processes maintain a simple form of the ion 
velocity distribution, 
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from which we calculate the corresponding fluorescence 
spectrum and fit the measurements, using a and w as 
free parameters. L(S±,w) oc (w 2 + S 2 .)^ 1 represents the 
Lorentzian lincshape of the optical pumping transition, 
with S± = kv z ± App. Equation ^ follows from the 
steady state of the underlying optical Bloch equations, 
augmented by a simple Krook-type collision term [3] , and 
smoothly interpolates between the collisionless (a = 0) 
and the strong-collision (a>l) regime [36] . As shown 
in Fig. [2^,, our data is well described by the fitted velocity 
distributions, f^\v z ) (Fig. [^d), from which we extract 
the average ion velocity v z = j v z f+ t \v z )dv z (Fig. [2^). 

We have also performed quantum-classical simulations 
of the laser-driven plasma dynamics and subsequent re- 
laxation in order to confirm our analysis procedure and 
to extend the parameter range of our study. The calcu- 
lations track the ion motion via classical molecular dy- 
namics (MD) simulations of the plasma ions, interacting 
by bare Coulomb interactions in a cubic simulation cell 
with periodic boundary conditions. The internal states 
of the ions are propagated according to the optical Bloch 
equations [37] corresponding to the laser-driven 4-level 
scheme shown in FigfTja. This allows to follow the time 
evolution of the internal-state density matrix during opti- 
cal pumping alongside the phase-space trajectory of each 
individual ion. The quantum spin evolution is coupled 
to the ion velocities via the Doppler shift, such that this 
approach yields the fully correlated internal and transla- 
tional plasma dynamics. The good agreement with our 
measurements (Fig. [2J:) demonstrates that this approach 
captures the essential physics of our experiments. In ad- 
dition, we used the simulated distributions to generate 
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FIG. 3: (color online), (a) Simulation results for the time 
evolution of the average ion velocity v z (symbols), obtained 
for a narrow perturbation of the initial velocity distribution 
centered around = v th = \Jkv,T/m. Equations Q and 
fitted to the data in the domain t < 3u> p t (lines), provide 
an excellent description of the short time dynamics. Panel (b) 
shows the extracted relaxation rate as a function of v z . At 
small F it displays a strong velocity dependence, which, how- 
ever, weakens dramatically in the strongly coupled regime. 
The black line shows a comparison to the weak-coupling 
Landau-Spitzer prediction (see equation |l])) for F = 0.005. 



theoretical spectra, which are then analyzed just as the 
experimental spectra are. The excellent agreement be- 
tween both results (Fig{2j;) attests to the accuracy of our 
experimental approach for extracting the average ion ve- 
locity based on equation 

To determine the underlying relaxation rate, we model 
the relaxation dynamics by a non-Markovian damping 
term EE] 



dt 



MM 0) ,t')v z (t-t')dt' , (4) 



where the stationary memory kernel M z (v z °\ t) accounts 
for retardation effects due to the strongly coupled na- 
ture of the equilibrium plasma. Here, ±ui°^ are the 
velocities on resonance with the pumping lasers, which 
show the strongest deviation from a Maxwellian. The 
deviation is assumed to be well-localized in this deriva- 
tion. The memory time vanishes in weakly coupled 
plasmas, and M z (v Zl t) = 2j z (v z )5(t). This yields fa- 
miliar exponential relaxation with a damping constant 
lz{v z ) = J dv±j(v) / eq (v_i_) that coincides with the 
Landau-Spitzer result, averaged over the corresponding 
Maxwellian of the transverse velocity vj^ = (v x ,v y ). 

However, temporal correlations become important for 
the strong coupling conditions of our experiments. Fol- 
lowing |39j . a simple Gaussian memory kernel 
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properly accounts for short-time correlations. The cor- 
responding memory time r is connected to the average 
ion acceleration [40] and can be obtained independently 
from equilibrium simulations. For the time scales rele- 
vant to our experiments (cj p t < 3), this simple theory is 
well confirmed by our simulations (Fig. [3]). 

An additional complication may arise from the veloc- 
ity dependence of the relaxation rate (cf. equation Q), 
which can make the dynamics of v z depend on the specific 
form of /+(v z ). To investigate this point we have per- 
formed simulations where the initial state transfer is done 
within a much narrower velocity range (Av z = 0.05 Vth) 
centered around a velocity ±v z °\ Fitting eqs.Q and 
|s| to our simulation results for different vi ^ yields the 
velocity-dependent rates shown in Fig. [3}d. The veloc- 
ity dependence weakens dramatically with increasing T 
and nearly vanishes for our experimental conditions. In 
fact, the optical pumping affects almost the entire ex- 
tent of the initially Gaussian velocity distribution (Fig. 
2}d). Since j(v) varies by less than 10% over this range, 
we can apply equations Q and ([s]) to our measurements 
and identify the extracted relaxation rate as the aver- 
age 7 ps J 7(w)/ eq (v)dv. As Fig. [4] demonstrates, this 
simple approach provides an excellent description of our 
measurements . 

Figure [5] summarizes our main results, showing the re- 
laxation rate 7 for a wide range of parameters. Upon 
scaling the rate by the ionic plasma frequency w p and ex- 
pressing the plasma temperature and density in terms of 
r, all data collapse onto a single universal curve. In the 
weakly coupled regime, the numerical results approach 
the familiar Landau-Spitzer form 
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7 = ar2 Wp lnA, A = 



(6) 



of the relaxation rate, with a = 0.46 and b = 0.53. Well 
into the strongly coupled regime, where equation ^ pre- 
dicts negative relaxation rates, we find good agreement 
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FIG. 4: (color online). Measured time evolution of the aver- 
age ion velocity v z for two different sets of plasma parameters, 
(a) p = 3.1 • 10 9 cm- 3 , T = 3.3K and (b) p = 4.4 • 10 9 cm -3 , 
T = 1.7K, corresponding to T = 1.2 and V = 2.6, respectively. 
The lines show a fit to the non-Markovian relaxation model, 
equations Q and |H]), from which we extract the correspond- 
ing relaxation rate 7 shown in Fig. [5] 
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FIG. 5: (color online). Average relaxation rate as a function 
of coupling strength I\ The large circles show the experimen- 
tal results for different combinations of plasma density and 
temperature, within a range of 6x 10 8 cm~ 3 < p < 5x 10 9 cm -3 
and 0.8K < T < 2.8K. Using dimensionless quantities (7/<^ P 
and F), the data collapses onto a universal curve, which veri- 
fies the expected Coulomb scaling and provides experimental 
evidence that, in the strongly coupled regime, the velocity 
dependence of the relaxation rate is negligible (see Fig|3]D) 
within our measurement accuracy. The thick solid line is ob- 
tained from MD simulations for a wider range of Coulomb 
coupling parameters. In the weak-coupling limit it approaches 
the Landau-Spitzer form, equation (pj, shown by the thick 
dotted line. The other lines show different proposed exten- 
sions [91412] into the strongly coupled regime, obtained by 
replacing A by A in equation (JsJ) according to the expressions 

given in the figure. The function E\(x) = f^^-j—dt denotes 
the exponential integral. 



between experimental and numerical results. A slight in- 
crease with r is evident in the experimental results and 
is well reproduced by our calculations. As T increases, 
both approach the plasma frequency, which sets the typ- 
ical time scale for ionic motion. 

Figure [5] also includes recently proposed theoretical ex- 
pressions based on effective Coulomb logarithms that ex- 
tend the Landau-Spitzer formula into the strongly cou- 
pled regime. We note, however, that the scaled rate 
will depend on the mass ratio of the species considered, 
and previous theory has focused on ion-electron thermal- 
ization with an eye towards dense plasma applications. 
Although our experimental accuracy appears sufficient 
to discriminate between different models, the additional 
mass-dependence, presently, limits such comparisons to 
a qualitative level. Perturbative corrections due to finite 
mass ratios have been investigated recently [H]. Exten- 
sions of existing models to equal-mass systems can now 
be subject to stringent tests through measurements in 
ultracold neutral plasmas. 

The described pump-probe technique makes a whole 
new class of experiments possible. Laser heating and 
cooling [35] will greatly stretch the range of accessible 
Coulomb coupling parameters and allow exploration of 
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the transition from an ideal to a correlated plasma, ex- 
tending more deeply into strongly coupled fluid regime. 
With improved time resolution, our approach will pro- 
vide experimental access to velocity autocorrelations and 
self-diffusion coefficients |43j . which determine dynamic 
structure factors and various transport processes. 
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